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Dedicated to Leonhard Euler on his 296th birthday 
1. Introduction 

Sometime around 1740, Euler || f|, |5j took a decisive step to unraveling the functional 
equation of the Riemann zeta function when he discovered a marvelous method of calculating 
the values of the (absolutely!) divergent series 



"1 + 1 + 1 + 1 + 1 + • 
"1 + 2 + 3 + 4 + 5 + - 



1 

~2' 
1 

~12' 

"1 + 4 + 9 + 16 + 25 + = 0, 

"1 + 8 + 27 + 64 + 125 + ■ • • " = — , etc.. 

In modern terms, these are the values at non-positive integer arguments of the Riemann zeta 
function ((s), which is defined by the series, absolutely convergent in Re(s) > 1: 

j- / \ -i 1 1 1 1 
v 1 2 s 3 s 4 s 5 s 

Naturally, having no notion of the analytic continuation at that time, to say nothing of 
functions of complex variable, Euler had to find a way of giving a meaning to those values 
of divergent series. What he actually did proceeds as follows. First, he directs his attention 
to "less divergent" alternating series 

r» _ 2 m + 3 m - 4 m + 5 m - 6 m + 7 m - 8 m + etc. 

since its convergent counterpart 

11111111 

D 1 1 1 h etc. 

2 n 3 n 4 n 5 n Q n 7" 8 n 

does indeed have faster convergence and is linked with the original series by the simple 
relation 

(i) as) = (i-2^x( S ) 1 



Work in part supported by Grant-in Aid for Scientific Research (B) No. 11440010, and by Grant-in Aid 
for Exploratory Research No. 13874004, Japan Society for the Promotion of Science. 



1 



where we have put 

~ — — -1-1 

*W ~ 2 s 3 s 4 s 5 s 

Then, he observes that the value © is obtained as a "limit" of the power series 

(2) l m - 2 m x + 3 m x 2 - 4Tx 3 + 5"V - 6 m x 5 + etc. 

as x — > 1, since, although the series itself converges only for |x| < 1, it has an expression as 
a rational function (analytic continuation, as we now put it), finite at x — 1, which is ob- 
tained by a successive application of multiplication by x and differentiation (or equivalently, 
applying the Euler operator x ■ d/dx successively after once multiplied by x) to the geometric 
series expansion 

(3) — — = 1 - x + x 2 - x 3 + x 4 - x 5 H , (\x\ < 1). 

1 + x 

For instance, if we substitute x — 1 in @, we find formally 

i = i-i + i-i-i-... = C(o) 

and hence, in view of ([[]), we have ((0) = —1/2. A few more examples are 

2x + 3x 2 - 4x 3 + 5x 4 , 

2 2 x + 3 V - 4V + 5 V , 

2 3 x + 3 s x 2 - 4 s x 3 + 5V , etc.. 

C(-2)=0, f(-3) = ~--- 

C(-2) = 0, C(-3) = 4-. 

For all that splendid idea however, this method unfortunately provides no rigorous way 
to establish the values of ((— m) as values of the analytically continued function ((s) at 
s = —m. (This may only be an afterthought and merely shows the degree to which Euler 
was ahead of his time and how much our modern point of view owes to him.) 

In the present article, aiming to evaluate the value £(— m) in as elementary, and yet 
rigorous, way as possible as the value of analytically continued function C( s ) ; we introduce 
and investigate a new g-analogue of the Riemann zeta function. As becomes clear in the 
course of our study, this function serves very well for the purpose not only of computing 
£(— m) but also of giving a nice g-analogue of ((s) valid for all s e C. 



which give us 



and in turn 



1 

(1 + x) 2 

1 — X 
{1 + x) 3 
1 — Ax + x 2 
(1 + x) 4 

c(-i; 

C(-i) = 



i 

4' 
1 

12" 
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To be more specific, as an alternative for a series like (H), we put x = g* (note x ■ d/dx 
is essentially d/dt) and, instead of repeating differentiation (this inevitably restricts us to 
looking only at the integer arguments), we replace n m by the g-integer [n] q := (f — g n )/(f — g) 
raised by the power — s (recall that Euler is the grand "Master of g"!); namely, we consider 
the series 



Throughout the paper, we always assume < q < 1, so the series (|) converges absolutely 
for any s 6 C and Re(i) > 0. If Re(s) > 1 (and Re(t) > 0), the series obviously converges 
to £(s) when q f 1. This suggests that we should regard the function f q (s, t) as a g-analogue 
of the Riemann zeta function ((s), but we reserve this until we make the specialization 
t — s — 1 which turns out to be utterly crucial. Before going into the specialization, we 
establish below the meromorphic continuation of f q (s, t) as a function of two variables s and 
t, which is carried out quite easily by using the binomial theorem. 

In the next section, we specialize t — s — 1 and establish a formula for s = —m G Z< 
(Proposition 2) as well as its limit when q f 1 (Theorem 1). Then we give the result 
concerning the limit as q f 1 for any s (Theorem 2). 

Proposition 1. Let < q < 1. As a function of (s,t) G C 2 , f q (s,t) is continued meromor- 
phically via the series expansion 

having poles of order 1 att & Z<o + 2iriZ/ logg := {a + 2nib/ logg | a, 6 G Z, a < 0}. 

Proof. We just apply the binomial expansion (1 — q n )~ s = ( s+r_ )g Tir and change the 

order of summations to get 

n=l v ^ ' n=l r=0 ^ ' 

= (i-g)'E( a + r " 1 )f;g B( ^ ) = (i-g)'f:( a + r " A ^ 



r=0 x 7 n=l r=0 

The other assertions follows readily from this 



r I 1 — g i+r 



Remark. It is worth noting that the function f q (s,t) can be expressed as the (beta-like) 
Jackson integral. In fact, we have 

q- t (l-q) 1 - s f q (s,t) = (l-q)J2 h Ws = / ^(l - gx)"%x. 

J=0 1 1 q ) Jo 
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2. Main results 



Now we put £ — s — 1. When s = —m G Z< , the point (s,t) = (— m, — m — 1) lies on 
the pole divisor £ = — m — 1 of f q (s,t). Nevertheless, a sort of "miracle" happens that the 
point turns out to be what is called "the point of indeterminacy", the function f q (s, s — 1) 
having a finite limit as s -> — m and moreover the limit approaches to the "correct" value 
£(— m) as g | 1. What is more, the function f q (s, s — 1) converges as g j 1 to for an?/ s ! 
These results, to be proved in quite elementary ways (certainly with only devices of which 
Euler could avail himself), well justify the function f q (s, s — 1) being referred to as the "true" 
g-analogue of the Riemann zeta function, and we label it hereafter as 



X) q n(s-l) q s-l q 2(s-l) q 3(s-l) q 4(s-l) 



M s q [i]; [2]; [3]; [4]- 

Remark. 1) Proper choice of £ seems to be essential. For example, the choice t = s adopted 
in m needed an extra term to adjust the convergence when q | 1 and gave no nice values at 
negative integers. The choices £ = s — 2, s — 3, s — 4, . . . seem as good as the value £(— m) is 
concerned, but extra poles at s — 2,3,4, ... emerge. However, these poles disappear at the 
limit q j 1. For example, with £ = s — 2 the residue at the simple pole s = 2is— (1 — q) 2 / logg 
which goes to as q | 1- How things become different depending on the choice of £ still 
seems to be mysterious. 

2) If we introduce the g-analogue ( q (s) of the alternating ((s) in the introduction by 

7 n(s— 1) 



COO = D" 1 



the identity corresponding to ([!]) takes the form 

C q (s)=( q (s)-2(l + qy s C q2 (s). 

In contrast to the situation of Euler, this does not help much and indeed even makes things 
worse because of the occurrence of another base g 2 . It may be said that once g is introduced, 
the acceleration of convergence is fully achieved and nothing more is needed. 

The formula in Proposition 1 when specialized to £ = s — 1 becomes 



(5) cm = d-d-EfT 1 ) q 



s+r— 1 



1 - g^- 1 

g*' 1 g s s(s + 1) g s+1 

q) 1 1 g-i + 1-g* 2 l-gs+i + 



Proposition 2. 1) The function ( q (s) has a simple pole at points in 1 + 27uZ/logg and in 
the set {a + 2nib/\ogq \a, b G Z, a < 0, b ^ 0}. In particular, s = 1 is a simple pole of ( q (s) 
with residue (g — l)/logg. 
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2) Form G Z, m > 0, the limiting value lim s ^_ m ( q (s) exists (which we write ( q (—m)) and 
is given explicitly by 



(6) 



C„(-m) = (!-„)- {£(-!)' (7) f 



\m+l 



(m + 1) logg 



Proof. Assertion 1) is straightforward from the formula lim 2/ _ > o|//(l — g y ) = — 1/logg 
being used for the residue at s = 1. For 2), note the terms with r > m + 2 in the sum 
vanishes when s — > —m since ( _m \j/* _1 ) = and 1 — q- m + r - 1 0. On the other hand, for 
r = m + 1 we have lim s ^_ m (s + m)/(l — g s+m ) = — \j logg and hence 

lim ( s + m \ qS+m (-l) m m\ / 1 \ (-1)^+1 
s^-m\m+ljl — q s+m (m + 1)! \ logg/ (m+l)logg 

The rest of the computation is clear. I 

Before giving our general formula for lim^i C?( — m ) (with expected value), let us look 
at the first few examples. 

Example 1. As stated in Proposition 2, ( q (s) has a simple pole at s — 1 with residue 
(g — 1)/ logg, which converges to 1 as q — > 1. This agrees with the well-known fact (reviewed 
later) that ((s) has a simple pole at s = 1 with residue 1. 

Example 2. By (0) we have 

1 1 

C*(o) 



q — 1 log g 
Since 

+ - + 0(g-i; 



logg log(l + (g-l)) (g-l)-(g-l) 2 /2 + --- g-1 2 
we find 

limC,(0) = -i. 

q— >i z 

This agrees with Euler's computation £(0) = —1/2. 
Example 3. Again by @ we have 

C,(-l) = (1 - 9)" 1 f ^-r - ^-r + ) 

\ q 2 — 1 g — 1 2 log q J 

1/1 1 1 1 



1 — q \q— 1 2 + g — 1 g — 1 2 log g 
1/1 lg-1 1 1 1 g-1 



l-gV.2(g-l) 4 8 g-1 2(g-l) 4 24 

as g — > 1, 



12 

in accordance with Q(— 1) = —1/12 
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Let the Bernoulli numbers Bk be defined by the generating series 



te 



X t k 



k=0 



kY 



First several values are 

Bq = 1, B 1 = -, B 2 — -, B 3 = 0, B 4 = - — , B 5 = 0, B e = — , _B 7 = 0, . . . 
Now the general formula is the following 
Theorem 1. For each non-negative integer m, we have 

lim ( q (—m) = 



Bm+l 



<in m + l 

Proof. On account of formula (@), we have to show 



,,„ ;(1 _ 9) - m {g ( _ ir ( : )_L_ 



+ 



-1)' 



m+l 



B 



m+l 



(m + 1) logg | m + l 

(Note here that since the sum on the left is finite so we may replace the limit q f 1 by 
q — > 1.) Multiplying both sides by (— l) m+1 (m + 1) and changing r — > m + 1 — r, we see this 
is equivalent to 



-l) m 5, 



m+l- 



Writing 
and using 

we have 



1 1 r log g 



<f - 1 r e rl °s9-l logg 



T = £(-im 



k=0 



kV 



m 



m 



r=l x ' 

m+l / \ oo 

+ 1 >£(- 1 >V! 1 );£(- 1 >^ 

r=l V y k=0 



[r\ogq) k 1 
£;! logg 



r k {~l) k B k 



oo / m+l 

= E E(-u 

fc=0 \r=l 

Since the inner sum on the right can be calculated as 



m + l 
r 



(logg) 



fc-i 



k\ 



m+l / 
r=l ^ 



m + l 

r 



x^) ((i-xr +i -i) 



-i 



if fc = 0, 

if < k < m + 1, 

(_l)m+i( m + 1)! if A; = m + 1, 
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we find 

m+l 



[m 



+ l)£(-lW m = + B m+l (\ogqT + 0((loggr +1 ) (as g - 1). 

^ Vr - V q r - 1 log g 



From this and the expansion logg = q — 1 + 0((q — l) 2 ) (g — > 1), we obtain the desired 
result. I 

Remark. In view of Theorem 1, we may define the q- Bernoulli number B m (q) by 

B m {q) := -m( q (l -m) (m > 1). 

By (P) (letting m — > m — 1 and r — >• m — r) we have the closed formula 

*•<*> - ta- 1 )--{|:(-ir(:)^T+^ 
= (s-d^'D-d'Pt^t (m - 1) ' 

Here, the term with r = in the last sum should be read as 1/ logg (the limiting value when 
r — > 0). This suggests to put 

logg 

With this, the g-Bernoulli numbers {-Bm(g)} m>0 satisfy the recursion 



m=0 v 7 



g m £? m (g) = (-l) n 5 n (g) + 5 ln {n > 0), 



where 5\ n = 1 if n = 1 and otherwise, and the generating function 

t m 



F q (t) :=J2 B m(Q) ^ 

m=0 

satisfies the relation 

F q (qt) = e l F q (t) - te l . 

This g-Bernoulli number is essentially (i.e., up to the sign (— l) m ) the same as the one 
introduced in Tsumura @. 

The following fundamental relation, apart from its own importance, guarantees that our 
computation at negative integers above does give us the correct values which we intended to 
obtain on a rigorous basis. 

Theorem 2. For any s G C ; s ^ 1, we have 

\im( q (s)=as)- 

<?Ti 
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What we understand as the right-hand side for arbitrary s is the value of the function 
analytically continued to the whole s-plane. (We give the analytic continuation by using the 
Euler-Maclaurin summation formula, see the proof below.) On the left-hand side, q should 
avoid the values with which ( q (s) has a pole at s, but this is achieved once q gets close 
enough to 1. 

Example. We give some numerical examples. Take s — 1/2 and q = 0.999 in (^). Sum of 
the first 10 5 terms gives us the value —1.46014527395 ■ • • . Replacing q by q = 0.99999 and 
taking the first 10 7 terms we get the value —1.460352417- • •, which agrees with the actual 
value C(l/2) = -1.4603545088 • • • up to 5 decimal points. Take the point s = 1/2 + 14.1347i 
near the first non-trivial zero (= 1/2 + 14.134725141734693790457251983562 •• -i) of ((s). 

For q = 0.9999, the first 10 5 terms gives the absurdly large 10835.552 h 10270.785 • • -i, 

while 10 6 terms gives -0.000306477 ■• • + 0.000794677 ■■■ i (the actual value is ((1/2 + 

14.1347z) = 0.000003135364 0.00001969336 ••■ i). If we take s = 1/2 + 14.134725i 

and q = 0.99999, the first 2 ■ 10 6 terms gives -0.4690527 0.4669811 •• -i and the 

5 ■ 10 6 terms -0.000031064 • • ■ + 0.0000812513 • • • % (the actual value is C(l/2 + 14.134725z) = 
0.000000017674 0.00000011102 ■■■%). 

Combining Theorem 1 and Theorem 2, we readily obtain 

Corollary. For each non-negative integer m, we have 

B m +i 



<r" + 2 m + 3 m + 4 m + 5 m + • • • " = C(-m) 



m + 1 



Remarks. 1) We can also define a g-analogue of the Hurwitz zeta function £(s; a) 

C,(s;a) = ^ — 

and prove the identity 



^ q (n+a)(s-l) 

\n + a\„ 



limC g (s;a) = ((s; a) 
for any s ^ 1, as well as the formula 

.B m+ i(a) 



limC,(-m;a) = 
?Ti m + 1 

for integers m < 0. Here, the Bernoulli polynomial B^ix) is defined by the generating series 

fc=0 

As in the remark after Theorem 1, we can also define the g-Bernoulli polynomial similarly 
and derive elementary formulas. But to make our presentation as concise as possible, we 
restrict ourselves to the case of the Riemann zeta function. 
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2) It would be amusing to note that the limit 



oc 



(8) Hm(l-g)*5]^^ = (A ; -l)!C(*) (Vfc>2,fc6Z) 

is derived easily from 



n=l 



limC g (fc) = C(fc). 



(The latter directly follows from the definition without appealing to Theorem 2 because we 
are in the region of absolute convergence.) In fact, if we put s = 2 in ([|) and make r + 1 — > n, 
we have 



C 9 (2) = (l-g) 2 ^ T 



nq 



n 



n=l H 



which gives the desired limit for k = 2. For general k, we similarly put s = k in (]5|) and 
make fc + r — 1 — > n to find 



&(*) = (!-*)*£ U-i 

n=l ^ ' 

(Observe ( fc "J = for n = 1, 2, . . . ,k — 2.) We note that 

+ lower degree terms 



h — 1 



Jc-lJ {k-l)\ 

and, on taking the limit q j 1, sums coming from lower terms vanish inductively, hence we 
obtain the conclusion. 

When is even and k > 4, the series X/nLi ^ fe ~ 1 5 ,n /(l — q n ) = Y^n=i (^2d\n^ k ~j 9™ 
constitutes the Fourier series of the Eisenstein series Gfc(r) of weight k on the modular 
group, with constant term — Bj s /2k(= ((1 — k)/2). Here r is a variable in the upper- 
half plane and is linked with q by q = e 2mT . The modularity amounts to the transformation 
formula G^—l/r) = r Gfc(T), which can be derived from, as Hecke showed, the functional 
equation of the corresponding Dirichlet series (p(s) := ((s)((s + 1 — k): 

(2tt)-T(sMs) = (-l) k/2 (2Tt) s - k T(k - s)<p(k - s). 

(When k is odd, the functional equation of (p(s) fails to take this form and so the series 
Y^Li n h ~ 1 q n /(l — q n ) cannot be a Fourier series of a modular form.) Hecke also showed that 
the residue of ip(s) at the simple pole s = k is equal to (2ni) k co/(k — 1)! where c is the 
constant term of the corresponding modular form. In our case, the residue is ((k) and thus 
the constant term of Gkij) is {k — l)\£(k)/( y 27ii) k = —Bk/2k, as expected. As an alternative 
way, we may use (^) to determine the constant term as follows: Put r = it with t > 0. Then 
e 2*i(-i/*t) _^ o as t -> and so 

1 (it) k 
c = JimG fc (--) = Vm(it) k G k {it) = lim-^(l - q) k G k (it) 

= — -^(k- l)K(k). 
(2vri) fcV ; sv ; 
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Proof of Theorem 2. Recall the celebrated summation formula of Euler |], |[ (and Maclaurin, 
cf. m §7.21], obtained simply by repeating integration by parts): For a C^-function f(x) 
on [1, oo) and arbitrary integers M > 0, N > 1, we have 



N „at M 



/ 1 
/(*) rfx + -(/(l) + /(AO) + £ (/W(JV) - 

»=1 fc=l ^ 

- WTTy.l B ««W/ <M+1) (-)^. 

where Sj^ + i(x) is the "periodic Bernoulli polynomial" defined by 

Bk{x) = Bk{x — [x]) ([x] is the largest integer not exceeds x). 

Recall the Bernoulli polynomial B^{x) is defined by the generating series (0): 

1 1 3 1 

B (x) = 1, Bi(x) = x - -, B 2 {x) = x 2 - x + -, B 3 (x) = x 3 - -x 2 + -x, . . . 

2 o 2 2 

As is well-known, by taking f(x) = x~ s and letting A" — ► oo, we obtain the analytic contin- 
uation of £(s) to the region Re(s) > — M: 



M 



(10 ) cw ^4 +E J^ Wi -i£k±y - Bm+i(x)x 



k ^(k+\y" (M+i)!^ 



-s-M-1 



where (s)^ := s(s + 1) • • • (s + k — 1). Since we may choose M arbitrary large, this gives 
the analytic continuation of £(s) to the whole s-plane, revealing the (unique) simple pole at 
s = 1 with residue 1. 

Now we take f(x) = q x ^ s ~^/(l — q x ) s and M — 1 in (|9|). Assuming Re(s) > 1 and 
noting 

fix) = \oz< r <r s ! ' 



and in general f^ k \x) = (\ogq) k q x( - s ~ 1 \l — q x y s ~ k x (a polynomial in s and q x ), we see that 
we can take the limit A" — > oo and obtain 

00 qn(s— l) poo qx(s—i) -y q s ~ l \ s — 1 + q 

5(13^ = / (r^ rf -+2-(T^-^ (log,), *"(r^ 



(1 ° M)2 r ft(j) ^ «(' + !) -3.(1 + 

(1 — g I ) s + 2 



for Re(s) > 1. The first integral on the right is computed as 

qX (s-i) d ^_r co q~ x _ Uq~ x - If-'V* 3 _ q s - l {l-q) l - s 
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x n(s-l) 



We therefore obtain 

00 

(11) C Q (s) = (i-^Ett^ 

n=l ^ ' 

g 8-1 q - 1 g 5 ' 1 g 5 ' 1 log? , l , 

s-1 ' logg 2 12 g-r ?j 

_ (1 _ qY (J2Ml j™~ BMf ^ ± ± 1) - y i - y (i - <r? dx . 

Unlike the classical case flTO|), just to let M larger does not make the convergence of the 
integral better, since the factor g x(s_1 ) in f( M+l \x) always forces Re(s) > 1. Instead, we 
use in (|Tl~|) the Fourier expansion of the periodic Bernoulli polynomials^ (cf. ||, Ch.IX, 
Misc. Ex. 12]) 

p 2mnx 

(12) B k (x) = -k\ 

The equality is valid for all real numbers x when k > 2, the sum being absolutely and 
uniformly convergent. Putting this (for k = 2) into (|IT ) and interchanging the summation 
and the integration, we find 

s — 1 logg 2 12 q — 1 

Wl ,2 V- 1 f° :E f s -i)S(s + l)-3s(l-g :E ) + (l-g*) 2 , 

+ 1-g logg 2 > - , / e 2 ™V (s 1} - ,/ x ' — dx. 

^ (2mn) 2 J x (1 - q x ) s + 2 

Further we make a change of variable q x = u to obtain 

< 13 » «.) - £L.tl + Cl + ^(.-i + „ 



1 - <?) s logg X] r9 1 a 2 i s ( s + l ) h ^ s -l + 8n,-s- 1) 



(2nin 

ngZ v 
n#0 



— 3s6 g (s — 1 + 8n, —s) + b q (s — 1 + 5n, — s + 1)} 
where we have put 5 = 2ni/ logg and 

bt{a,p) = [ u a -\l -uf~ x du 



(referred to as the incomplete beta functionfj]). Note that each of the incomplete beta integrals 
in flT3| ) converges absolutely for Re(s) > 1 and uniformly bounded with respect to n; 



a — 2/-r „,\ — (J+v — 1, 



\b q (s - 1 + 8n, -s + v)\ < / m^^I-u)^ 1 ^ (Vn, <7 = Re(s), v = -1,0,1) 



*We owe Ueno-Nishizawa || the idea of replacing £2 0*0 in the integral by its Fourier expansion. However, 
our argument that follows, which uses only integration by parts and no confluent hypergeometric functions 
or the like, seems considerably different from the one in ||. 

^We are tempted to remind the reader that the beta integral is often called the Euler integral. 
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hence the sum converges absolutely. 

Now, repeated use of integration by parts provides us the formula 

lh{a,P) = f (— \ (l-uf- l du = -t a (l-t)> 3 - 1 -1—^- f u a {l-uf- 2 du 
Jo \ a J a a Jo 

= U a {l-tf- 1 -^—^- f (^—\ (l- u f- 2 du 
a V ; a J \a + lj K ! 



\k-l (1-/9) k— 1 i_a+k— 1 



J2(-i) k - i y ^ k - l r +k -\i-tf- k 



k=i 



>M-1 (! _ P)m-i 



+(-l) M ~ l v . C Pt{a + M-1,P-M + 1) 
{a)M-i 

for any M > 2. Applying this to b q (s — 1 + 5n, — s — 1), we have (note q Sn = 1) 

M-i , , 

b q (s-l + 5n,-s-l) = ^(-1)^-1 [S _ +Z l k -\ ^(l-g)— 1 -* 

+ (£ + ^_ 6 (s — 2 + M + <5n, — s — M). 

(s-l + <Sn) M -i 

This accomplishes the analytic continuation of 6 g (s — 1 + 5n, — s — 1) as a function of s into 
the region Re(s) > 2 — M. From this we have 



(27rm 

n^O 



M-l , s 

^ K ' ^ (2nin) 2 (s - 1 + 8n) k q 1 q> 



k=l ™ez 



+ (_i) M -'y-L- — ^k±i — [ 9 u s - 3+M+Sn (i- u ) 



-s-Af-1 



M-l 



^ V - (s)fc+1 q s+k - 2 (l - q)- s - 1 - k 

— L> (2mn) 2 (s - 1 + 5n) k q { q) 



k=l 



- (-l^-Mogg^ 1 e ^n XqX{ s-2 + M ){l _ qX) -s-M-l 

•f^f (27rm) 2 (s - 1 + dnjM-i A 

Using 



lim^- = -l, lim(l-g) fc (s-l + 5n)fc = (-2«) fc , lim ^—^ = ^ 

1 — q q-+\ g->l 1 — q 
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we obtain, for Re(s) > 2 — M, 



s(s + 1) 

lim(l - q) s logg V , . , 2 b q (s -l + 6n, -s - 1) 



nez \ 



^ (27rm) fc + 2 ^ (27rm) M+1 A 

k=l ngZ v / ngZ v I J 1 



M-l 

^(* + 2)l' 

In the last equality, we have used flT2| ) and its specialization {x = 1) 



(27rm) fc fc! 

n#0 



nSZ 



valid for all k >2. We do exactly the same for the terms containing b q (s — 1 + Sn, —s) and 
b q (s — 1 + 5n, — s + 1). As it turns out however, the contributions from these two vanish when 
we take q | 1, for the powers of 1 — q involved are lower than those from b q (s — 1 + 5n, — s — 1) . 
We therefore obtain, for Re(s) > 2 — M, 



1 , 1 , B k+ i (s)m+1 f°° u ( \ -s-M-l j 

= h > t; rr [s )k — ~, rr / B M+1 [x)x ax. 

s-l 2 ^{k + iy} ! (M + 1)\J 1 + v ; 

This coincides with formula flTTi| ) for ((s) valid in Re(s) > —M, and thus the theorem is 
established since the integer M can be arbitrary large. ■ 
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